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We present a comprehensive review of available crater topography measurements for Saturn’s moon 
Titan. In general, the depths of Titan’s craters are within the range of depths observed for similarly sized 
fresh craters on Ganymede, but several hundreds of meters shallower than Ganymede’s average depth vs. 
diameter trend. Depth-to-diameter ratios are between 0.0012 ± 0.0003 (for the largest crater studied, 
Menrva, D ~ 425 km) and 0.017 ± 0.004 (for the smallest crater studied, Ksa, D ~ 39 km). When we eval- 
uate the Anderson-Darling goodness-of-ht parameter, we find that there is less than a 10% probability 
that Titan's craters have a current depth distribution that is consistent with the depth distribution of 
fresh craters on Ganymede. There is, however, a much higher probability that the relative depths are uni- 
formly distributed between 0 (fresh) and 1 (completely infilled). This distribution is consistent with an 
infilling process that is relatively constant with time, such as aeolian deposition. Assuming that Gany- 
mede represents a close ‘airless’ analogue to Titan, the difference in depths represents the first quantita- 
tive measure of the amount of modification that has shaped Titan’s surface, the only body in the outer 
Solar System with extensive surface-atmosphere exchange. 

© 2012 Elsevier Inc. All rights reserved. 


1. Introduction 

Unique among the icy satellites, Titan’s surface shows evidence 
for extensive modification by fluvial and aeolian processes 
(Tomasko et al„ 2005; Lorenz et al., 2006, 2008; Stofan et ah, 
2007). These processes act to change the topography of its surface 
over time, through fluvial erosion, mass wasting, burial by dunes 
and submergence in seas. Quantifying the extent of this modifica- 
tion is difficult, since the original, un-eroded surface topography is 
generally unknown. However, fresh craters on icy satellites have 
well-known shapes and morphologies, which have been deter- 
mined from extensive studies of the airless worlds of the outer 
Solar System (e.g., Schenk et ah, 2004). By comparing the topogra- 
phy of craters on Titan to similarly sized, relatively pristine ana- 
logues on airless bodies, we can obtain one of the few direct 
measures of the amount of modification that has occurred on Titan. 

The best analogues for comparison to Titan are Jupiter’s moons 
Ganymede and Callisto. For large craters (formed in the gravity re- 
gime), crater size is dependent on gravity, impact velocity, projec- 
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tile size, and target and projectile density (Holsapple and Housen, 
2007). Ganymede and Titan have similar gravity (g ~ 1.4 m/s 2 ), tar- 
get density (p ~ 1 g/cm 3 ), and likely projectile density (p — 1 g/ 
cm 3 ), and the average impact velocity at Ganymede is twice that 
at Titan (201<m/s vs. 10.5 km/s, Zahnle et ah, 2003). The gravity 
on Callisto is also similar to that of Titan (g= 1.25 m/s 2 ), and its 
average impact velocity is somewhat lower than that on Gany- 
mede «v) - 15 km/s). However, given that the depths of craters 
on Callisto are almost indistinguishable from those of Ganymede 
(Schenk, 2002), and the topographic data set that exists for Callisto 
is somewhat sparser than that for Ganymede, we focus in this pa- 
per solely on comparisons to Ganymede. 

In addition to target and projectile properties, the subsurface 
structure of an icy satellite is important to the shape and morphol- 
ogy of craters, since rheological changes at depth influence the 
depth-to-diameter ratio of large craters (D > 26 km) (Schenk, 
2002). This range of diameters coincides with those most likely 
to be observed on Titan, as models predict that craters with 
D > 20 km were formed by projectiles large enough to be only 
minimally disrupted by Titan’s extended, thick atmosphere 
(Korycansky and Zahnle, 2005). In the case of Ganymede and Titan, 
it appears that their outermost layers (~1001<m) are both 
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dominated by cold water ice. Titan’s long-wavelength topography 
(Nirnmo and Bills, 2010) and Cassini Radio Science Subsystem 
(RSS) gravity measurements (less et al., 2010) both suggest that 
Titan has a floating, isostatically compensated ice shell with a 
mean thickness of ~100km overlying a subsurface ocean. 
Ganymede is also thought to possess an H 2 0-NH 3 ocean at 
~100 km depth (Spohn and Schubert, 2003). 

Cassini RADAR has imaged ~50% of the surface of Titan at reso- 
lutions as good as 350 m (e.g., Elachi et al., 2005; Hayes et al., 

2011) , and in this data set, more than 60 potential craters have 
been identified (Wood et al., 2010; Neish and Lorenz, 2012). Topo- 
graphic information for these craters has thus far been difficult to 
obtain. There is only a limited amount of topographic data on Titan 
in the form of altimetry (Zebker et al., 2009) and stereo (~2% areal 
coverage, Kirk et al., 2012). Additionally, altimetry requires nadir 
pointing, and thus is incompatible with simultaneous synthetic 
aperture radar (SAR) imaging of the surface. Fortunately, a tech- 
nique known as ‘SARTopo’ (Stiles et al., 2009) has been developed 
that estimates surface heights by comparing the calibration of 
overlapping SAR beams. The Cassini RADAR instrument has five 
different antenna feeds (or beams), and data from overlapping 
beams are acquired nearly simultaneously. Surface height can be 
computed by maximizing the correlation between the received 
power and the antenna gain pattern at each point along track in 
the overlap region. This technique is capable of estimating surface 
heights for most of the SAR-imaged surface of Titan with ~10 km 
horizontal resolution and a vertical resolution of tens of meters. 
It has extended the area over which co-located topography and 
SAR imagery is available on Titan by at least an order of magnitude. 

In this work, we present topography data for several craters on 
Titan using the SARTopo data set. We compare this topography to 
similarly sized craters on Ganymede, for which topography has 
been extracted from shadow length measurements and stereo- 
derived digital elevation models (Schenk, 2002; Bray et al., 2008, 

2012) . Finally, we make inferences regarding the relative amount 
of landscape degradation that has occurred on Titan due to erosion 
and infill. 

2. Observations 

We plotted the position of the SARTopo data over SAR images of 
every known crater on Titan with D > 20 km. There are roughly 30 
such craters (Wood et al., 2010; Neish and Lorenz, 2012), with 
eight interpreted as ‘certain’ impact structures. Of these, six ‘cer- 
tain’ or ‘nearly certain’ craters (Table 1) and two ‘probable’ craters 
had corresponding SARTopo coverage (Fig. 1). Of the six ‘certain’ 
craters, five had topographic profiles broadly consistent with 


craters of their size on other icy worlds, including a central peak 
in the 39 km diameter crater Ksa (Kirk et al., 2012; see also 
Figs, lb and 2), a central uplift in the 100 km diameter crater Hano 
(Figs, lg and 3), and a near-flat profile for the 425 km diameter cra- 
ter Menrva (Figs, li and 3). Only Soi lacked any recognizable topog- 
raphy (Fig. 4), as did the two ‘probable’ craters (craters #43 and 
#49 from Wood et al. (2010); Fig. 5). Both of the ‘probable’ craters 
are located in Titan’s expansive sand seas, so these features are 
either not impact craters or have been completely infilled with 
aeolian deposits. Since we cannot unambiguously determine 
whether these features are craters, we exclude them from the 
remainder of the depth analysis. 

The diameter of each crater was determined from the SAR imag- 
ery, and the errors reported in Table 1 represent the natural varia- 
tions of the craters away from circularity. Height measurements 
were determined from a 10-km long (along-track) by w-km wide 
(across-track) region of SAR pixels, where w is the width of overlap 
between the two antenna beams that contribute to a SARTopo pro- 
file. The value of w varies depending upon the spacecraft altitude 
and instrument pointing angles, but is generally on the order of 
several kilometers. All of the pixels in the 10-by-w km region were 
then used together with knowledge of the viewing geometry and 
radar antenna gain patterns to estimate the heights. The along- 
track translation between consecutive measurements is signifi- 
cantly smaller than the 1 0 km region of interest (typically 
200 m), however, making it necessary to subsample tracks to a 
10 km separation distance (~50 data points) to ensure indepen- 
dent height measurements. This was possible for all depth 
measurements reported here. Note that this method of height 
derivation tends to reduce the calculated crater depths, since 
crater rims are averaged with adjacent, lower topography and cra- 
ter floors are averaged with adjacent, higher topography (see 
Fig. 6). These measurements could therefore be considered lower 
limits. 

For the six ‘certain’ craters, we then calculated depth, 
d = hi - h 2 , by taking the difference between the highest point on 
the crater rim and the lowest point on the crater floor, on both 
sides of the crater, d i and d 2 (Fig. 6). Systematic errors in height, 
Shj, were propagated throughout the analysis. These errors were 
determined from radar instrument noise and viewing geometry, 
and are summarized in Table 1 of Stiles et al. (2009). For the depth 
measurements, we used the topographic profile closest to the mid- 
point of the crater. However, in several cases, only profiles near the 
edge of the crater were available (e.g., Sinlap, Afekan). The depth 
measurements therefore assume that the depth is constant across 
the crater floor. If this assumption is incorrect, these values repre- 
sent lower limits for the depth of the craters. 


Table 1 

Depth-to-diameter ratio for seven ‘certain’ or ‘nearly certain’ craters on Titan. 


Crater 

Diameter, D (km) 

Depth, d (m) 

d/D 

Technique 

Relative depth, R b,c 

Relative depth, R e 

Ksa 

Momoy 

Soi 

Sinlap 

Hano 

Afekan 

Menrva 

39 ±2 

40 ±1 
78 ±2 
82 ±2 

100 ±5 
115 ± 5 
425 ± 25 

660!™ 
750 ±175 
680 ±100 
110±100 
640^ 
700 ±100 
525+Jf 
455 

490 +12° 

0 01 7+ 0 004 

U.U 1 /_o.o 04 
0.019 ±0.005 
0.017 ±0.003 
0.001 ± 0.001 
°008 US® 

0.009 ± 0.001 
o nn5 +0001 

U.UUD_q qoi 

0 OO4 +0 002 
U.UU‘i_o 002 

0 001 ?+°’° 003 
u.uuiz_ 00003 

SARTopo 

Stereo 3 

Autostereo 

SARTopo 

SARTopo 

Autostereo 

SARTopo 

SARTopo 

SARTopo 

0 94+ 0 - 20 

U.Z^_ o.20 

0 1 3+ 0 - 20 
u 1J -0.20 

022+CJJ 

0.901SS 

0.43!™ 

0.38!“;“ 

0 5fi +0 ' 09 
U.DO_oo8 

0.62!“;]| d 

N/A 

0 49+ 015 
O.^Z_015 

0 ^4+ 015 
0.40!°;“ 
0.89!™ 

fl qc+0.16 

0-36_ 015 

0.30!“;]“ 

0.46!]];]; 

0 S7 +0 ’ 19 
U.DZ_ 01 g 

N/A 


a See Kirk et al. (2012). 

b Relative depth is defined as R(D ) = 1 - d t (D)/dg(D), where d t (D ) is the depth of a crater with diameter D on Titan, and dg(D) is the depth of a crater with diameter D on 
Ganymede. 

c Ganymede crater depths from Table 4 in Bray et al. (2012). 
d Assumed to have the same depth as a D = 100 km crater. 
e Ganymede crater depths from Fig. 2b in Schenk (2002). 
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Fig. 1 . Cassini RADAR images of craters on Titan with depth estimates: (a) ‘Probable’ crater #43 from Wood et al. (2010), D ~ 26 km. (b) ‘Certain’ crater Ksa, D ~ 39 km. (c) 
‘Certain’ crater Momoy, D ~ 40 km. This crater was discovered in Titan flyby T77, located near 11.6°N, 44.6°W. Note the putative impact melt flows just to the east of the 
crater, (d) ‘Probable’ crater #49 from Wood et al. (2010), D ~ 60 km. (e) ‘Certain’ crater Soi, D ~ 78 km. (f) ‘Certain’ crater Sinlap, D ~ 82 km. (g) ‘Nearly certain’ crater Hano, 
D ~ 100 km. (h) ‘Certain’ crater Afekan, D ~ 115 km. (i) ‘Certain’ crater Menrva, D ~ 425 km. Each image is presented in simple cylindrical projection, and north is up. 
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Fig. 2. (top) Ksa crater (D ~ 39 km) as viewed by Cassini RADAR. North is up, and illumination is from the bottom, (bottom) Three topographic profiles through Ksa crater. The 
light gray points represent the data and the black lines the uncertainty in their along-track and vertical position. The position of the profiles is shown in the RADAR image 
above. The size of each point is determined by its associated along-track and across-track errors. 
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Fig. 3. Topographic profiles through five craters with corresponding SARTopo data and topography consistent with craters on other icy satellites: Ksa (D ~ 39 km), Sinlap 
(D ~ 82 km), Hano (D ~ 100 km), Afekan (D ~ 1 15 km) and Menrva (D ~ 425 km). The top four craters are plotted at the same spatial scale covered by Menrva for comparison. 
Representative along-track and vertical error bars are included in each plot. 
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Fig. 4. (top) Soi (D ~ 78 km), as viewed by Cassini RADAR. North is up, and illumination is from the bottom, (bottom) Topographic profile through Soi. The light gray points 
represent the data and the black lines the uncertainty in their along-track and vertical position. Gray arrows highlight the position of the western and eastern rims of the 
crater apparent in the Cassini RADAR image. 
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Fig. 5. (top) ‘Probable’ crater #43 (D ~ 26 km) identified by Wood et al. (2010), as viewed by Cassini RADAR. North is up. and illumination is from the bottom, (bottom) 
Topography near the probable crater. The light gray points represent the data and the black lines the uncertainty in their along-track and vertical position. No topography 
consistent with a crater is visible in this region, suggesting that the feature is not an impact crater or that it has been completely filled by sand. 
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Fig. 6. Topographic profile of Ksa crater, marked with the location of the rim and 
floor elevations used to determine crater depth. The depths on either side of the 
crater, d and d 2 , were averaged to account for the natural variation of crater depth 
along the rim. Each point in this data set was averaged over a~10x2 km wide 
area. This method tends to reduce crater depth, as adjacent topography is averaged 
with the rim maximum and floor minimum. 


The two depths, d ^ and d 2 (where d ^ is taken to be the larger 
depth and d 2 the smaller) were then used to determine the average 
depth of the crater. We recognize that natural variations in depth 


around the rim of a crater, which are known to vary by ~10-20% 
on Mars (Stewart and Valiant, 2006), will likely dominate the 
uncertainty in the crater depth. Thus, to be as conservative as pos- 
sible in our depth estimation, we report a range based on the larg- 
est depth consistent with our measurements, d ^ + Sd^, and the 
smallest depth consistent with our measurements, d 2 - Sd 2 . Stereo 
coverage is necessary to fully quantify the variations in depth 
around the rim of the crater, but thus far, Cassini RADAR has only 
acquired a stereo pair over Ksa crater. Preliminary analysis of this 
pair gives a depth of 750 m, although the elevation of the rim var- 
ies by as much as 350 m, representing a variation in depth of ~45% 
(Kirk et al., 2012). 

In the case of Soi, which was deemed a ‘certain’ crater by Wood 
et al. (2010), the SARtopo profile is nearly flat (Fig. 4; the position 
of the crater rim apparent in the SAR image is denoted by arrows). 
From morphological evidence in the SAR image, we judge that Soi 
is an impact crater (albeit, a highly modified one), and assign a 
depth based on the position of the inferred crater rim. We find that 
the western edge of Soi has a depth of 220 ± 60 m, whereas the 
eastern edge is essentially flat. There is also a topographic high 
near the middle of the crater, possibly a central peak similar to that 
observed at Ksa (Kirk et al., 2012). Averaging the depths on the two 
sides of the crater, we assign Soi a depth of 110 ± 100 m. 

Depths were also estimated for two of Titan’s craters by com- 
paring the foreshortening of the near and far walls, in a technique 
we refer to as ‘autostereo’, which assumes a perfectly symmetric 
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profile across the crater. The depth of one of these craters (Sinlap, 
diameter D ~ 82 km) agrees well with the corresponding SARTopo 
data, supporting the idea that Sinlap has a relatively flat floor. (An 
earlier estimate of 1300 m reported in Elachi et al. (2005) was late- 
ly found to be in error by a factor of two, resulting from confusion 
about the pixel scale of the image measured.) The second crater 
(Momoy, D ~ 40 km) was only recently discovered in RADAR imag- 
ery of the surface of Titan from Cassini flyby T77 (Fig. lc). It is 
roughly the same depth (680 m vs. 660 m) as the similarly sized 
Ksa crater (D ~ 39 km). 


3. Comparisons to Ganymede 

Topography data were acquired for seven ‘certain’ or 'nearly 
certain’ craters on Titan (Table 1). The depth-to-diameter ratios 
for these craters span the range 0.0012 ± 0.0003 (for the largest 
crater studied, Menrva, D~425 km) to 0.017 ±0.004 (for the 
smallest crater studied, Ksa, D ~ 39 km). We can compare these 
depths to fresh craters on Ganymede, which we assume represents 
a close ‘airless’ analogue to Titan (with few exogenic processes 
modifying the surface). We find that the depths of Titan craters 
are generally within the range of depths observed for similarly- 
sized craters on Ganymede, but several hundreds of meters shal- 



Fig. 7. Depth of craters on Titan (gray circles) derived from SARtopo measurements 
compared to similarly sized, fresh craters on Ganymede (central peaks, squares; 
central pits, stars) (a) from Bray et al. (2012) and (b) Schenk (2002). Dashed lines 
represent the power law fits to the Ganymede data reported in those two papers. 
Error bars in the Bray et al. (2012) data set incorporate both the range of natural 
variation in crater depth and a systematic error of ~5%. As shadow length 
measurements are only made from one quadrant of each crater, the natural 
variation in crater depth cannot be determined for the Schenk (2002) dataset. We 
do, however, include systematic errors of 30%, which represent the maximum 
uncertainty in the shadow length measurements due to residual uncertainties in 
the control network and viewing geometry, as well as the limitation that the 
measurements can only be done on one side of the crater (Schenk, 1991). 


lower than the average depth vs. diameter trend for Ganymede 
(Fig. 7). 

Given these data, we used the Anderson-Darling goodness-of- 
fit technique to evaluate the likelihood that the crater depths on Ti- 
tan represent fresh Ganymede craters that have been altered by 
infilling. This technique evaluates the statistic A 2 , which quantifies 
the difference between the observed data and model distribution 
functions, F n (x) and F(x), respectively: 

A 2 = nj [(F„(x) - F(x)) 2 /(F(x)(\ - F(x)))]dx (1) 

The greater the difference between the data and model probability 
distribution, the larger A 2 becomes, increasing the probability of 
observing a given A 2 value by chance. In the case where the model 
F(x) has no free parameters to adjust, the probability, 
p = 1 - P(A 2 < z), of observing a particular A 2 value can be obtained 
from standard tables (Lewis, 1961 ). Here, P is the probability distri- 
bution function of the Anderson-Darling statistic described in Lewis 
(1961), and z is the independent variable in that probability distri- 
bution function. A high p value indicates that the observed result is 
likely if the null hypothesis being tested is true, whereas a low p va- 
lue indicates a significant contradiction to the null hypothesis. 

Using this method, we evaluated two different hypotheses: (1) 
that the Titan craters we observe are fresh and have similar depths 
to fresh Ganymede craters of the same size and (2) that the ob- 
served depth of Titan’s craters has been influenced by erosional in- 
fill after the crater was formed. We compared Titan depth 
measurements to Ganymede’s fresh crater depths from two differ- 
ent data sets: the depths of craters on Ganymede reported in Bray 
et al. (2012), calculated from Galileo stereo DEMs (Fig. 7a), and the 
depths reported in Schenk (2002), calculated primarily from sha- 
dow length measurements (Fig. 7b). The depths in Bray et al. 
(2012) are notably shallower than those calculated by Schenk 
(2002), plotting parallel to, but below, the trend line in Schenk 
(2002). This difference can be explained by the averaging of the 
stereo data, which causes a decrease in measured crater depths. 
Similarly, the averaging of the Titan topography data over a 
10 km x wkm region may reduce crater depths compared to the 
higher resolution Ganymede data (which have resolutions of sev- 
eral hundred meters to more than a kilometer per pixel). However, 
when we compare the depth determined from the Ksa DEM (Kirk 
et al., 2012) at resolutions similar to those on Ganymede 
(1.4 km/pixel) to the depth obtained from SARTopo (~10 km/pix- 
el), we find little difference (750 ± 175 m vs. 660 ±170 m). We thus 
judge that resolution alone cannot account for the several hundred 
meter difference in depth observed in the Ganymede and Titan 
data sets. 

We now evaluate the two hypotheses listed above. Hypothesis 1 
states that Titan’s craters were selected from the depth distribu- 
tion of fresh craters on Ganymede. For the Bray et al. (2012) data 
set, this hypothesis has a significance of A 2 = 2.03, or p = 0.09, if 
one only considers the upper limits of the crater depths for those 
craters with ‘recognizable’ topographic profiles. This significance 
drops to zero if we include Soi, or consider only the lower limits 
of the crater depths, since no known fresh crater on Ganymede 
with D >20 km has a depth less than ~500 m. For the Schenk 
(2002) data set, the significance is zero in either case. Hypothesis 
2 states that the relative depths of Titan’s craters are uniformly dis- 
tributed between 0 (fresh) and 1 (completely infilled). Here, the 
relative depth, R, is given by R(D) = 1 d t (D)/dg(D), where d t (D) is 
the depth of a crater with diameter D on Titan, and d g (D) is the 
depth of a crater with diameter D on Ganymede, calculated from 
the empirical relation for d determined by Schenk (2002) and the 
relation for d/D determined by Bray et al. (2012). This hypothesis 
has a much higher significance than Hypothesis 1, A 2 = 0.931, or 
p = 0.4, for the Schenk (2002) data set, and A 2 = 0.330, or p = 0.92, 
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for the Bray et al. (2012) data set. If we only consider lower limits 
on depth, this hypothesis has a significance of A 2 = 1.470, orp = 0.2, 
for the Schenk (2002) data set, and A 2 = 0.949, or p = 0.38, for the 
Bray et al. (2012) data set. If we only consider upper limits on 
depth, this hypothesis has a significance of A 2 = 1.341, or p = 0.22, 
for the Schenk (2002) data set, and A 2 = 0.851, or p = 0.45, for the 
Bray et al. (2012) data set. These significances are increased further 
if we assume a uniform distribution between 0 and x, where x < 1 
as a consequence of there being a minimum detectable depth of 
~100 m in the SARtopo data set. 


4. Discussion 

Given the results of our analysis, we judge that it is very unli- 
kely that Titan's craters have the same depth distribution as fresh 
craters on Ganymede. They are more consistent with a uniform 
distribution of relative depths between 0 (fresh) and 1 (completely 
infilled). Such a distribution suggests a modification process whose 
rate is constant with time, such as aeolian infilling (Forsberg- 
Taylor et al., 2004). This interpretation is consistent with the loca- 
tion of the craters studied in this work, which are all (except Hano) 
within 30° of Titan’s equator, a region dominated by large sand 
seas (Lorenz and Radebaugh, 2009). In contrast, fluvial modifica- 
tion has an infilling rate that diminishes with time, as crater slopes 
diminish and the area of deposition on the crater floor increases. 
Craters in the Sinus Sabeus Quadrangle on Mars show a depth 
distribution consistent with fluvial modification, with a peak near 
R = 0.05 (fresh craters) and another peak near R = 0.6-0.8 (de- 
graded craters) (Forsberg-Taylor et al„ 2004). In this case, R is 
defined as R = (H - h)IH, where H represents the original fresh cra- 
ter depth as estimated from relationships derived from MOLA 
topography (Bleacher et al., 2003) and h is the actual crater depth. 
Only one crater with R ~ 0.6-0.8 was observed in our data set, 
suggesting either that fluvial modification is not the primary 
erosive process on Titan, or that the craters we are observing are 
very young and have not yet reached the state of degradation 
observed on Mars. 

A uniform distribution of relative depths is also consistent with 
direct atmospheric sedimentation, but we judge this hypothesis to 
be less compelling given the expected photochemical production 
rates on Titan. The most abundant solid product of Titan’s photo- 
chemistry is predicted to be acetylene, C 2 H 2 , and photochemical 
models predict a downward flux of C 2 H 2 between 5.6 x 10 s and 
1.9 x 10 9 molecules cm~ 2 s _1 (Yung et al., 1984; Lara et al., 1994; 
Toublanc et al., 1995). This flux represents a sedimentation rate 
(assumed uniform over Titan) of ~10-35 m/Gyr. Even if Titan’s cra- 
ter retention age is older than the 1 Ga that current estimates sug- 
gest (Neish and Lorenz, 2012), it is unlikely that a layer of directly 
precipitated atmospheric sediment could account for the several 
hundred meters of infill we interpret to be present in Titan’s 
craters. 

We thus judge aeolian infill (which may represent horizontal 
surface transport of photochemical material that was originally 
deposited from the atmosphere) to be the more likely cause. This 
would represent an infilling rate of several hundred meters of sand 
per billion years, if sand transport rates have been constant on Ti- 
tan over these timescales. If the dunes represent a younger surface 
feature, the inferred infilling rate could be much higher. On Earth, 
the floor of the crater Roter Kamm appears to have been filled by 
more than 300 m of sand on timescales less than a few million 
years (Fudali, 1973; Radebaugh et al., 2010). Thus, many more cra- 
ters might be present but unidentified on Titan due to rapid aeolian 
infilling. 

Another mechanism known to alter crater topography and re- 
duce crater depths on icy satellites is viscous relaxation (Passey 


and Shoemaker, 1982; Dombard and McKinnon, 2006). Craters 
>10 km in diameter on Ganymede show a range of relaxation 
states, from fresh craters to craters with more subdued topography 
and upbowed floors. However, given Titan’s lower surface temper- 
ature (95 K vs. 120 K), viscous relaxation is predicted to cause less 
than a 3% change in topography for craters with D < 125 km on Ti- 
tan (Baugh and Brown, 2006), suggesting that viscous relaxation is 
not an important mechanism for modification of most craters on 
Titan. Nonetheless, craters near the polar regions of Ganymede, 
where surface temperatures approach 80 K, do show evidence for 
relaxation (Dombard and McKinnon, 2006). In addition, recent 
measurements of the depths of large (D > 50 km) craters on Rhea 
and Dione suggest that numerical models under-predict the 
amount of crater relaxation, indicating that these satellites might 
have experienced more heating early in their history than previ- 
ously thought (Phillips et al., 2012). 

Still, we judge it unlikely that viscous relaxation has been the 
primary form of crater modification on Titan. Viscous relaxation, 
like fluvial modification, has a modification rate that diminishes 
with time (see, for example, Fig. 9 in Dombard and McKinnon, 
2006). This is because the strain rate associated with viscous flow 
is proportional to the depth of the crater (Melosh, 1989). As the 
crater relaxes and the depth decreases, the rate of modification will 
slow. We would therefore not expect a uniform distribution of rel- 
ative depth values if viscous relaxation were the primary mecha- 
nism for crater modification on Titan. In addition, we would 
expect smaller craters to have systematically smaller relative 
depths, R, than larger craters, since larger craters relax faster than 
smaller craters (Passey and Shoemaker, 1982). This relationship is 
not observed for the Schenk (2002) data set, but further work will 
be required to determine whether it is consistent with the Bray 
et al. (2012) data set (Fig. 8). Of the craters with D < 125 km stud- 
ied in this work, Afekan seems most consistent with a relaxed cra- 
ter, showing evidence for an upbowed crater floor in the crater 
profile (Fig. 3). 

Unlike its smaller counterparts, the largest crater observed on 
Titan, Menrva (D ~ 425 km), would be expected to undergo signif- 
icant topographic relaxation (assuming it did not breach the crust 
altogether, or that compositional differences such as the presence 
of methane clathrate do not make Titan’s near-surface rheology 
substantially different from Ganymede’s). On Ganymede, craters 
with D>~150km are significantly shallower than smaller 
complex craters (Schenk et al., 2004). The largest preserved impact 
feature on Ganymede, Gilgamesh (D ~ 590 km), shares a similar 
topographic expression to that of Menrva. Gilgamesh has a central, 
domical surface rising up to ~500 m and spanning ~150km, 



Fig. 8. Distribution of relative depth with crater diameter, calculated relative to the 
Ganymede depths from Bray et al. (2012) (gray) and Schenk (2002) (black). 
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similar to that observed at Menrva (Fig. 3), though with an addi- 
tional ~1 km high scarp bounding the dome (Schenk et al., 
1997). Gilgamesh’s dome is then surrounded by a 225 km wide 
annular zone of hummocky terrain, bounded by an inward-facing 
scarp ~l-km high, which is presumed to be the basin rim. Beyond 
the apparent rim lies a second scarp 0.5-1 km high at a diameter of 
750 km, and a furrow at a diameter of ~980 km which is concentric 
to Gilgamesh. This morphology gives Gilgamesh a resemblance to 
multi-ring basins on the Moon, such as Orientale (Passey and Shoe- 
maker, 1982). Menrva’s rim is somewhat more subdued than that 
of Gilgamesh, with a height of only several hundred meters, per- 
haps eroded by the fluvial processes evident in this region. It also 
lacks the exterior scarps observed at Gilgamesh, and thus cannot 
be classified as a multi-ring basin. It is possible that if such scarps 
existed in the plains beyond Menrva’s rim, they would be buried by 
sediment and indistinguishable to Cassini RADAR. Numerical mod- 
eling of the formation and evolution of Menrva should be a produc- 
tive direction for future research. 

The Visual and Infrared Mapping Spectrometer (VIMS) provides 
compositional information for the crater floors on Titan, providing 
an independent means to distinguish between the different modi- 
fication processes. VIMS has observed three craters on Titan in de- 
tail: the 82 km diameter Sinlap (Le Mouelic et al., 2008), the 90 km 
diameter Selk (Soderblom et al., 2010), and the 120 km diameter 
Paxsi (Buratti et al., 2012). Spectra of the floors of both Sinlap 
and Selk indicate they are composed of the ‘dark blue’ spectral unit, 
which has been interpreted as being enhanced in water-ice (Rodri- 
guez et al., 2006; Barnes et al., 2007). Spectra of the floor of Paxsi 
indicate that it is composed of the ‘dark brown’ spectral unit, 
which is correlated with regions of dune fields (Soderblom et al., 
2007). Thus, only one of the three craters studied is consistent with 
the presence of aeolian infill (although only one of these - Sinlap - 
has the topographic data necessary for our depth analysis). The 
spectral properties suggest that Selk and Sinlap are fresh craters, 
or that they have been altered by viscous relaxation or fluvial mod- 
ification, rather than aeolian infilling. Fluvial modification is most 
consistent with their morphologies (see, for example, Fig. 5 in Fors- 
berg-Taylor et al. (2004)), and would act to deposit icy alluvium in 
the crater floor, consistent with their spectral properties. Taken to- 
gether, Titan’s craters seem more consistent with aeolian infill, but 
individual craters such as Sinlap may have experienced some flu- 
vial modification. 

In addition to the six certain craters discussed above, two ‘prob- 
able’ craters were found to have overlapping SARtopo data. Neither 
was observed to have topography consistent with a crater, but both 
features are located in Titan’s sand seas, leading to the possibility 
that they have been completely infilled by aeolian deposits. Simple 
aeolian infilling tends to leave the crater rim largely clear of depos- 
its, while the crater center becomes filled with a deposit roughly 
parabolic in cross-section (Forsberg-Taylor et al., 2004). Thus, we 
could observe a relatively pristine crater rim if these craters had 
been buried in sand. However, the topography data appear to be 
more consistent with a small hill around which sand has been di- 
verted. Future stereo measurements of one or more of the probable 
craters would aid in determining whether they were caused by an 
impact. 

Although we have assumed in this work that Titan’s craters rep- 
resent modified versions of Ganymede’s fresh craters, there are a 
few possible differences between the crusts of Titan and Ganymede 
that might affect crater morphology. Two of these are (a) the pres- 
ence of clathrates, such as methane clathrate (i.e., Tobie et al., 
2006), or (b) the presence of liquid hydrocarbons in the crust (i.e., 
Kossacki and Lorenz, 1996). Models of Titan’s interior evolution 
suggest that several of Titan’s volatile species may be stored in 
the crust as clathrates (Tobie et al., 2012). Indeed, the outgassing 
of methane clathrates in Titan’s upper crust has been suggested 


as a possible mechanism for replenishing Titan’s photochemically 
unstable atmosphere (Tobie et al., 2006). Experiments have shown 
that methane clathrate hydrate is very rigid compared to water ice. 
At temperatures just below the freezing point of water, methane 
hydrate is over 20 times stronger than water ice at a given strain 
rate, and the contrast increases at lower temperatures (Durham 
et al., 2003). In general, for projectiles with the same characteristics, 
stronger targets will yield smaller craters than weaker targets (e.g., 
Wiinnemann et al., 2011). So if Titan’s crust consists of a thick layer 
of methane clathrate, we may expect craters on Titan to be smaller 
than those formed on Ganymede by comparable projectiles. How- 
ever, models indicate that the clathrate layer, if present, should 
be relatively thin (<10 km) compared to the thick, water ice crust 
below (>50 km) (Tobie et al., 2006), although the details of enclath- 
ration and retention of porosity (which will also influence strength) 
have not been explored. This scenario - impacts into a strong layer 
of methane clathrate over a weaker layer of water ice - has not been 
studied with impact hydrocodes. However, simulations into lay- 
ered basalt show that the inclusion of thin layers of stronger basalt 
(two orders of magnitude stronger than normal basalt) does not 
have any significant effects on crater dimensions (Senft and Stew- 
art, 2007). If clathrates are not present on Titan, or are only present 
as thin layers, this should have a relatively minor effect on crater 
dimensions compared to those on Ganymede. Thick layers of clath- 
rates, however, may have an observable effect on crater diameter, 
and we thus judge this to be an area worthy of further study. 

The presence of liquids on the surface and in the near subsurface 
can also cause extensive modification to crater shape, as is observed 
on Earth (e.g., Collins and Wiinnemann, 2005; Ormo et al., 2006). In 
the case of Titan, liquids consist of hydrocarbons, either as wet 
unconsolidated sediments (such as those observed at the Huygens 
landing site; Niemann et al., 2005) or shallow marine environments 
(such as the lakes observed at the north and south poles; Stofan 
et al., 2007). Craters formed in similar environments on the Earth 
lack any significant surface topography, including the absence of 
a raised rim, as poorly consolidated, water-saturated sediments 
slump into the crater, generating a much larger structure than is ex- 
pected from the amount of kinetic energy released from the impact 
(Collins and Wiinnemann, 2005). In addition, sedimentation subse- 
quent to the impact has covered or subdued the topographic 
expressions of all known marine-target craters on Earth (although 
buried craters can later be exposed at the surface by erosion, i.e., 
Lockne crater in Sweden, Lindstrom et al., 1996). It is possible that 
the lack of topography associated with marine-target impacts may 
help to explain the relative dearth of impact craters observed near 
the poles of Titan (Wood et al., 2010). If Titan’s polar regions are sat- 
urated by liquid hydrocarbons, craters formed in those regions may 
lack any recognizable topographic expression. 


5. Summary 

In summary, our new observations of craters on Titan suggest 
that: 

(1) The depths of Titan’s craters are within the range of depths 
observed for similarly-sized fresh craters on Ganymede, 
but several hundreds of meters shallower than Ganymede’s 
average depth vs. diameter trend. 

(2) There is less than a 10% probability that Titan’s craters have 
a depth distribution that is consistent with the depth distri- 
bution of fresh craters on Ganymede. 

(3) The differences in depth appear most consistent with aeo- 
lian infilling, although some amount of viscous relaxation 
and fluvial modification may have contributed to the modi- 
fication of Titan’s craters. 
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